DFTT 06/2002 
ROM2F/2002/02 



Effect of the galactic halo modeling on the DAMA/Nal annual 
modulation result: an extended analysis of the data for WIMPs 
with a purely spin-independent coupling 

P. Belli, 1 '!] R. Ceruffi/'fl N. Fornengo, 2 '0 and S. Scopel 2 '| 

1 Dipartimento di Fisica, Universitd di Roma "Tor Vergata" and INFN, 
Sezione di Roma2, 1-00133 Roma, Italy 
2 Dipartimento di Fisica Teorica, Universitd di Torino and INFN, 
Sezione di Torino, Via P. Giuria 1, 1-10125 Torino, Italy 
(Dated: February 1, 2008) 

Abstract 



The DAMA/Nal Collaboration has observed a 4-er C.L. model independent effect investigating 
the annual modulation signature in the counting rate of an Nal(Tl) set-up (total exposure of 57986 
kg day) and the implications of this effect have been studied under different model-dependent as- 
sumptions. In this paper we extend one of the previous analyses, the case of a WIMP with a 
purely spin-independent coupling, by discussing in detail the implications on the results of the 
uncertainties on the dark matter galactic velocity distribution. We study in a systematic way pos- 
sible departures from the isothermal sphere model, which is the parameterization usually adopted 
to describe the halo. We specifically consider modifications arising from various matter density 
profiles, effects due to anisotropies of the velocity dispersion tensor and rotation of the galactic 
halo. The hypothesis of WIMP annual modulation, already favoured in the previous analysis using 
an isothermal sphere, is confirmed in all the investigated scenarios, and the effects of the different 
halo models on the determination of the allowed maximum-likelihood region in the WIMP mass 
and WIMP-nucleon cross-section are derived and discussed. 
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I. INTRODUCTION 



In the last few years various technical approaches have been exploited in direct searches 
for Weakly Interacting Massive Particles (WIMP) 0. In particular, the DAMA/Nal Col- 
laboration has collected a very large body of statistics which allows one to look for the 
distinctive signature in Dark Matter (DM) direct detection represented by the annual mod- 
ulation of the rate @, |J , an effect which is due to the rotation of the Earth around the Sun 
||. The analysis of the DAMA/Nal data after 4 years of running 0, |J, corresponding to a 
total exposure of 57986 kg day, has indeed led to the observation of an annual-modulation 
effect, which does not appear to be related to any possible source of systematics ||. This 
exciting result has been analyzed under different hypotheses on the properties of WIMP dark 
matter: purely spin independent coupling M, mixed spin-coherent interaction ||, inelastic 
dark matter [7| . In the case of purely spin-independent interactions, the annual modulation 
result has been shown to be compatible with a galactic halo composed, at least partially, by 
relic neutralinos in different classes of supersymmetric models || 

In the present analysis we consider the case of a WIMP with coherent interactions domi- 
nant over the spin-dependent ones. The analysis of the counting rate of any direct detection 
experiment may be done in terms of the WIMP mass my/ and of the quantity i&^\^° n \ 
where cr s ^ 1 ^ 1 r eon ' ) is the WIMP-nucleon cross section for scalar interaction and £ 1 is the frac- 
tional amount of local non-baryonic DM density which is ascribed to the WIMP responsible 
for the effect (£ < 1) pj, Performing a maximum-likelihood analysis of the data, the 
DAMA/Nal Collaboration has derived a region in the plane rriy/-^ &^^ on ^ which is com- 
patible at 3-cr C.L. to the observed effect of annual modulation. The properties of this 
region are sensitive to astrophysical inputs 0, [Uj, some of which deserve a reanalysis and a 



deeper insight. 

One of the main ingredients for the calculation of the expected rates is the distribution 
function of WIMPs in their six-dimensional phase space: F(f, v) d 3 r d 3 v (where the position 
vector r = (x,y,z) and the velocity vector v = (v x ,v y ,v z ) are defined in the rest frame of 
the Galaxy). Direct detection rates Rdet depend on the distribution function (DF) at the 
Earth position in the Galaxy: 

/(0)=F(i2o,u), (1) 

where Ro = (Rq, 0, 0) is the location of the Earth at a distance Rq ~ 8.5 kpc from the 
galactic center and along the galactic plane. It is therefore clear that an accurate calculation 
of the expected detection rates requires a knowledge of the phase-space distribution function 
F(f, v) d 3 r d 3 v. 

From the observational side, the most relevant piece of information coming from astro- 



1 In Rcf. § the same symbol indicates a different quantity: £ = pw/(0.3 GeVcm 3 ). 
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physics is related to the rotational velocity of objects bounded to the Galaxy: 

where G is the Newton's constant and M tot {r) denotes the total mass contained inside the 
radius r = \r\: 

M tot (r) = [ d 3 r'p tot (r'). (3) 

J r'<r 

The rotational velocity therefore depends on the total matter density distribution in space: 

ptot (r) = Pdm if) + pvis (r) , (4) 

where the DM density distribution function is: 

p DM (f) = J d 3 vF(r,v), (5) 

and p V is{f) represents the contribution to the matter density due to components other than 
the DM, like the disk and the bulge. The local values for the rotational velocity and for the 
DM matter density are denoted by Vq = v rot (Ro) and po = Pdm(Ro) and they represent two 
key parameters in the calculation of WIMP direct detection rates, as it will be discussed in 
the following. 

In order to calculate the DF of Eq. ([!]) one must invert Eq. (|5|) taking into account 
observational data. This problem is affected by degeneracies that cannot be solved without 
adding some piece of information. This explains why the velocity distribution represents 
one of the main source of uncertainty in the calculation of direct detection signals. 

The usual approach to this problem consists in assuming that the system has some 
symmetry and that the distribution F depends on the phase space parameters only through 
some integrals of motion (energy, angular momentum): this last condition automatically 
implies stationarity and that the Jeans' equations are verified [11]. The velocity ellipsoid 



cr jj =< V{Vj > may than be calculated as a function of the derivatives of the potential (which 
are related to the rotational velocity) by making use of Euler's equation [11]]. Physically, this 



corresponds to impose hydrostatic equilibrium between pressure and gravitational attraction. 

The most common and widely used example of such a procedure, and by far the simplest, 
is the isothermal sphere model. It consists in a spherical infinite system with a flat rotational 
curve, which automatically implies p(r) oc r~ 2 and the potential \& oc ln(r 2 ) . The DF may be 
easily worked out, and turns out to be a Maxwellian: f{y) oc exp(— 3v 2 / (2v 2 ms )) , where v 
i ' and v rms denotes the root mean squared velocity of the WIMPs. The isothermal sphere 
describes a self-gravitating gas of collisionless particles in thermal equilibrium, representing 
the highest entropy rearrangement of WIMPs in their phase space. A strong argument in 
favour of this last property is the "violent relaxation" model of Lynden-Bell which 



indicates that the violently changing gravitational field of the newly formed Galaxy may 
have lead the non-interacting WIMPs to thermal equilibrium. Hydrostatic equilibrium 
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and the assumption that the velocity ellipsoid is isotropic allows to calculate v rms through 
the relation: v^ ms = 3/2 v^ ot (R ). Due to its simplicity, the isothermal sphere model has 
become the "standard" assumption in the evaluation of DM expected rates, and has been 
used extensively in the literature, including the analysis of the DAMA/Nal modulation data 
|| . However many of its underlying assumptions (sphericity of the halo, absence of rotation, 
isotropy of the dispersion tensor, flatness of the rotational curve) are not strongly constrained 
by astrophysical observations. Moreover the isothermal sphere is strictly unphysical and can 
only represent the behaviour of the inner part of physical systems, since it has a total infinite 
mass and needs some cut-off at large radii. 

In light of the latest experimental data on WIMP direct searches, the issue of possible 
departures from the isothermal sphere model has gained interest and prompted several dis- 
cussions H, [K| |i~3fl . In the present paper we intend to analyze this issue in a systematic 
way, by employing a comprehensive set of self-consistent galactic halo models. Each model 
introduces a different degree of deviation from the simple isothermal sphere. We consider 
modifications in the velocity distribution function which are originated from a change of the 



gravitational potential or a change of the DM density profile |TJ, [L5|, [If], [L7], [18], [19], |20 
We classify the different models depending on the symmetry properties of the galactic halo: 
spherical potential/density profile with an isotropic velocity dispersion; spherical poten- 
tial/density profile with a non-isotropic velocity dispersion; axisymmetric models; triaxial 
models. For the axisymmetric models we also consider the possibility of having a co-rotating 
or counter-rotating halo. 

The plan of the paper is as follows. In section II we introduce the formalism for the 
calculation of direct detection rates and the annual modulation signal, and summarize the 
procedure used by the DAMA/Nal Collaboration to determine the annual modulation region 
in the plane mw-^&lc&ilv ^ f° r a P ure ly spin-independent interacting WIMP. In Section III 
we describe the halo models that we intend to discuss and introduce a naming scheme that 
will be used throughout the paper. Section IV is devoted to the discussion of the constraints 
on the dark halo of our Galaxy coming from available observational data. In Section V the 
annual modulation region is calculated in a systematic way for all the models previously 
introduced, and the results are discussed. Finally, Section VI is devoted to our conclusions. 



II. DIRECT DETECTION RATES AND ANNUAL MODULATION EFFECT 

The expected differential event rate of a WIMP direct search experiment is given, for a 
monoatomic detector, by the expression: 

where Nt is the number of the target nuclei per unit of mass, mw is the WIMP mass, w 
and f(w) denote the WIMP velocity and DF in the Earth frame (w = \w\), da/dE^ is the 
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WIMP-nucleus differential cross section and Er is the nuclear recoil energy. Notice that 
the detection rate is directly proportional to the local DM density p . The generalization of 
Eq.(||) to a diatomic detector, like Nal, is straightforward. 

The differential cross section is, in general, a sum of a coherent and a spin-dependent 
contributions. In this paper we consider only the case of a WIMP whose interactions are 
dominated by the coherent part. In this case the rate may be expressed in terms of the 



WIMP-nucleon scalar cross section, cr^^° a> , as: 



(nucleon) 

CIO. 

2 



da ( _da_\ K( e r) ( 1 + m w /m p \ a ( nuc i e0 n) (?] 

dE R ^{dE R ) coherent ^ E T » \l + m w /m N ) ^ ' U 

where m p and m^ are the proton and nucleus mass, A is the nuclear mass number, E R ax is 
the maximal recoil energy and E n (E R ) is the nuclear form factor for coherent interactions 



which may be parametrized with the usual Helm expression [pi 
The relation between the velocities v and w is given by: 

v = v e + w, (8) 

V(B =V Q + V(Broti (9) 

where v@ and v & denote the velocities of the Earth and the Sun in the Galactic rest frame 
(\v & \ ~ vq+12 km/sec), and v^ rot is the Earth's orbital velocity around the Sun (^©rotl =30 
km/sec). Projecting Eq.([|) in the galactic plane, one gets: 

l%l = l^ol + \v®rot\ cos7cos[u;(t - t )] (10) 

where 7 is the inclination of the plane of rotation with respect to the galactic one, u = 2%/T 
with T=l year, and to — 2 nd June corresponds to the day when the Earth's velocity is at 
its maximum. 

The change of reference frame of Eqs. (|§1,|9|,|T0"|) introduces through the DF f(w) a time 
dependence in the expected rate Rdet- In order to exploit this time-dependence to extract 
the modulated signal from the measured counting rates, we follow the maximal-likelihood 
procedure of Ref . [Q , to which we refer for a detailed discussion. This procedure allows to 
determine the region in the plane (m w - £cr s \^ 1 r eon ' ) ) which is compatible with the modulation 
signal. A lower bound on mw at the value mw = 30 GeV is applied, to conform to the 
analysis of Ref. |2|, [| . 

The data we analyze in the present paper refer to the full set of data released so far 
by the DAMA/Nal Collaboration @ (DAMA/Nal 0-4), including also the upper limit on 
^scalar' 011 '' obtained by the same Collaboration, as discussed in Ref. 0. The same values of 
the quenching factors and of the cut efficiences as in Refs. |2|, (§] are used. We stress that also 
the determination of upper limits is affected by the choice of the WIMPs DF. This means 
that also when confronting upper limits one has to specify the galactic halo models which 
has been considered in the calculation. The formalism introduced in the present paper may 
in fact be used also to quantify the uncertainty in the determination of upper limits from 
direct detection experiments. 
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III. HALO MODELS 

Let us turn now to the discussion of the galactic halo models and of the techniques used to 
calculate the velocity distribution function f(v). The different models are classified according 
to the symmetry properties of their matter density profile (or gravitational potential) and 
of the velocity distribution function. We define four classes: 

(A) Spherically symmetric matter density with isotropic velocity dispersion; 

(B) Spherically symmetric matter density with non-isotropic velocity dispersion; 

(C) Axisymmetric models; 

(D) Triaxial models. 

All the models which we describe in this Section, and that will be used in the rest of the 
paper, are summarized in Table O, where we also introduce a naming scheme that will be 
of practical use in the discussion. 



A. Spherically symmetric matter density with isotropic velocity dispersion 

The first class of models is represented by those with a spherically symmetric matter 
density p{r) = pir) and isotropic velocity distribution f{v) = f{v). These two conditions 
imply that the phase-space DF depends on the space and velocity variables only through 
the energy, which is an integral of motion: F(f,v) = F(e), where e = \l/(r) — v 2 /2 is the 
relative energy (per unit mass) of the WIMP and \I/ is the relative potential, related to the 



total density p to t through the Poisson's equation [[Tl|] : 

V 2 ^ = -AnGptot. (11) 

Notice that ptot refers to all the matter components of the Galaxy, like the disk, the bulge 
or the halo, as written in Eq.(|J). 

Once the total potential \I/(r) is known, the WIMP DF F(e) may be worked out by 
inverting Eq.([5[). A change of variables from r to \l/ allows to cast Eq.([5|) as [ |TT|j : 



4vr £ F(e)yj2(V-e) de = p DM (V) , (12) 



By performing a Laplace inversion on Eq.(|12|), one obtains the Eddington formula |TT 



Fie) = 1 d f d f£M^l f S= i 13) 

In Eqs. (|T^,|T^) the normalization of which is defined through Poisson's equation up to an 
arbitrary constant, is fixed by requiring that \&(oo) = 0. The velocity distribution function 
f(v) which enters the calculation of direct detection rates is then obtained as in Eq.([I|). 
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Eq.(|13D shows that the dark matter DF F(r, v) depends not only on the halo DM density 
distribution Pdm, but also on the density distribution of all the other galactic components. 
This has to be the case, since the DM particles move inside the gravitational potential gen- 
erated by all the matter that makes up the Galaxy A complete and rigorous determination 
of the DM DF will therefore require to model not only the galactic halo, but also all the 
other (disk, bulge) components of the Galaxy (see, for instance Ref. pBfl). 

However, WIMP direct detection is directly sensitive only to local properties of the 
Galaxy, and in particular the relevant parameters in the calculation of the detection rate are 
the local rotational velocity v and the local DM matter density po- I n fact, Vq is directly 
related to the average WIMP kinetic energy, which is relevant in the scattering process with 
the nuclei of the detector, and the change of reference frame of Eqs. which is crucial in 

determining the amount of annual modulation in the detection rate. Instead, po is a sort of 
normalization factor for the direct detection rate Rdet- The dependence of Rdet essentially on 
local parameters implies that a detailed modeling of the inner (r <C Rq) part of the Galaxy, 
where the disk and bulge components are more relevant and in general dominant over the 
halo, is not crucial for our analysis. Moreover, the matter density of non-halo components 
at the local position in the Galaxy (r = Rq) are no longer dominant with respect to the halo 
matter density at the same position (see, for instance Ref. f23|). The bulge, in fact, can 
be described by using a spheroidal density distribution which gives a sizeable contribution 
inside the first kpc from the galactic center, and it is truncated at about r ~ 2 kpc. The 
disk has an exponential distribution which in most of the models dies away at about 4 kpc 
from the galactic center. We will therefore assume in the following that in the outer Galaxy 
the dominant contribution to the matter density is given by the halo: 

Pvis < Pdm for r > R . (14) 

The only basic information which is required from the non-halo components is their contri- 
bution to the local rotational velocity v : 

v 2 o = v 2 rot (Rq) = -|- [M vis + M halo ] , (15) 

where 

M vis iha i o) = / dV p vis {ha i ){r') . (16) 

Jr'<R 

A maximal halo occurs when M V i S <C Mhalo'- in this case almost all the local rotational 
velocity is supported by the halo and the local DM density po gets its maximal value p™ ax 
compatible with the given vq. The opposite situation occurs when M V i s assumes its maximal 
value compatible with observations: in this case, the local rotational velocity gets the max- 
imal contribution from the non-halo components and only a fraction of Vq is supported by 
the halo. At the same time po gets its minimal value p nm , for the same Vq. The constraints 
on these parameters are discussed in Section [TV]. 
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From the point of view of calculating the DM DF (Eq.fll3|)), the occurrence of a maximal 
or non-maximal halo modifies the gravitational potential ^f(r) and therefore the velocity 
distribution function f(v) is affected. Indicating with ^o(r) the potential for the maximal 
halo, the condition of Eq. (|14]) allows the generalization to the non-maximal 

¥(r) = -^-*o(r) + ( 1 - ^ vl . (17) 

Po \ Po / r 

The condition of Eq. ([14]) allows to work out the total potential for the case of non- 
maximal halos without explicitly modeling the visible parts of the Galaxy (bulk, disk): all 
the dependence of these components are contained in M V i S . 

As a comment, we notice that the presence of a non negligible contribution to vq from the 
non-halo components alters also the velocity distribution function of the isothermal sphere. 
The standard Maxwellian form for the isothermal sphere is in fact correct only for a maximal 
halo. 

Now that we have discussed the procedure to calculate the velocity distribution function 
once the matter density of the DM is given, we proceed to introduce the different models. 

The first type of model is a direct generalization of the isothermal sphere by introducing 
a core radius R c (model Al). The density profile is: 



2 ml + r 2 



PDM{r) - A^G (i£Vr>)2 ' (18) 
which corresponds to the following potential for a maximal halo: 

*o(r) = -f ln(i? 2 + r 2 ). (19) 

From the analytic form of this potential we will refer to this type of models as logarithmic 
models. The usual isothermal sphere (model AO) corresponds to the limit R c — > 0: 

which corresponds to the following potential for a maximal halo: 

vl/ (r) = -|ln(r 2 ). (21) 



For a maximal halo, the velocity distribution function arising from Eqs. (EG, 21 ) is the 



standard Maxwellian one. For a non-maximal halo, a deviation is present also in the R c — > 
as discussed above. The rotational curve supported by ^ of Eq. ( |19"D is: 

J ° (i?2 + r 2 



<t{r) = v 2 tt^—t, ■ (22) 



It is rising for small radii and at large radii becomes flat. In the Rc — > limit, v rot = Vq for 
all radii. 



S 



A second type of spherical models (A2 and A3) is defined by the following matter density 

P°m(t) - {Rl + r2){p+m ■ (23) 

which corresponds to the following potential for a maximal halo: 

Mr) = ^0yT 2 (^0). (24) 

From the analytic form of this potential we will refer to this type of models as power-law 
models. They represent the spherical limit of the more general class of axisymmetric "power- 
law" models of Ref. |16| which will be discussed in Section |IIIB| . The family of power-law 
models given by Eq.fl2"4]) is not defined for (3 = 0. However, when (3 = is substituted in Eq. 
(P3|), the density of Eq. fli"8|) is recovered. In fact the logarithmic model turns out to have 
the properties of the "missing" (3 = power law potential. Evaluating Eq. (|2~3"D for r = R 
the parameter \l/ a can be expressed in terms of the density po- The rotational velocity for 
the power-law models is given by: 
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2 _ y*g^ c r 
rot (_R2 + r 2U/3+2)/2 ' 

and it is asymptotically falling with r if (3 > (model A2) and rising if (3 < (model A3). 

The last family of spherical models we consider (models A4-A7) is defined by the following 
matter density: 

fR y\ l + (R /ar ] { ^ )/a 
P ™ = Po {t) [ 1 + (r/ar \ ■ (26) 

for the choice of values of the parameters a, (3, 7 and a summarized in Table |. Except 
for the Jaffe model (A4), the other three density profiles (A5, A6, A7) are obtained from 
numerical simulations of galaxy evolution. 



B. Spherically symmetric matter density with non isotropic velocity dispersion 

The procedure described in the previous Section can be generalized to the case of a non- 
isotropic velocity distribution, while keeping a spherically symmetric density profile. In this 
case, the most general DF is a function of e and of the magnitude of the angular momentum 
vector (L = \L\): F = F(e,L). Among the different choices which can be performed, a 
particularly simple case is obtained when F depends on e and L only through the so called 
Osipkov-Merrit variable |TT], EM : 

Q = "-^ w 

where the parameter r a is related to the degree of anisotropy (3q of the velocity dispersion 
tensor (evaluated at the Earth's position) in the following way 



7) 2 E>2 

A) = 1 - 4 = tfPhf- (28) 

V r Ro + r 2 a 
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Here the velocity is expressed in spherical coordinates and v </> = v g ^ v r (with vf =< vf > 
- < Vi > 2 , i = r,6,(p). 

The corresponding DF can be obtained by solving a modified version of Eddington's 
inversion formula, which is obtained by making the following substitutions in Eq . (|T3|) fll]: 



Q, (29) 
pDAiij) -> p Qt DAi(r) = ( 1 + ) Pdm{t) ■ (30) 



The models we consider are the same as discussed in Section [III A[ the logarithmic model 



of Eq. ( pT8|) (model Bl), the power-law models of Eq. ( pB"D (models B2 and B3) and the 
models defined by Eq.(^) (models B4-B7). The velocity distribution functions, which are 
obtained by solving the Eddington equation with the Osipkov-Merrit term, are therefore 
anisotropic with a degree of anisotropy controlled by the parameter /3 related to r a as in 
Eq.pi). 



C. Axisymmetric models 

In the case of axial symmetry the DF depends in general (at least) on two integrals of 
motion, the relative binding energy e and the component L z of the angular momentum along 
the axis of symmetry. The DF may be written quite generally as the sum of an even (F + ) 
and an odd (F_) contribution with respect to L z : 

F(e,L 2 ) = F + (e,L 2 ) + F_(e,L 2 ), (31) 

where: 

F ± = -[F(e,L g )±F{e,-L g )] . (32) 
When Eq.flT^p is extended to the axisymmetric case, the density pdm turns out to depend 



only on the even part F + jy], so that, by inverting it, the DF may be determined up to 
an arbitrary odd part F— The problem of the determination of F + for an axisymmetric 
matter density is both analitycally and numerically hard to perform and actually it requires 
a double Laplace inversion on Eq. (|T2"D . However, for particular families of axisymmetric 



potentials this problem has been solved analytically by Evans [15, 16 1. These families are 



the axisymmetric generalization of the first two classes already introduced in Section [HI A| : 



the first family has a logarithmic potential, the second one has a power-law potential. We 
stress that these analytic solution for F + are obtained under the assumption that the halo 
potential is dominant over the other components. Therefore these solutions correspond to a 
maximal halo. Even though they do not represent the most general situation, their simplicity 
makes them of practical use and convenient for studying the axisymmetric case. 
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The first family of axisymmetric potential we consider is the logarithmic potential fI5 
^models CI and C2): 



Vr\ , / „o „o Z 



%(R } z) = --±\n(Rl + R 2 +*-), (33) 
z q 

where R 2 = x 2 + y 2 is the radial coordinate along the galactic plane, R c is the core radius 
and q the flatness parameter. The corresponding DM density is: 

o (Rz)- vl (V + lR 2 + ^ 2 + (2-g-V (u) 

PDM(jR ' Z) ~ A^Gf (Rl + R* + z*q-*)* • (34) 

Eqs.(|33|) and ([34|) are the axisymmetric generalization of Eqs. ( |I9|) and (0). The corre- 
sponding rotational curve is obtained from Eq.(|22|) by substituting the radial coordinate r 
with the radial coordinate in the galactic plane R. 

By expressing z as a function of ^ through Eq . (|33|) , the density of Eq. (|34D may be 



decomposed as [16 



p = p (^) + R 2 Pl (^), (35) 
which allows to determine F + in the form: 

F + (e,L z )=F° + (e) + L 2 z Fl(e), (36) 

leading to a particularly simple analytic solution for the DF | I5| . We give it for completeness 
in Appendix |A|. The relation of Eq.(35) is no longer valid for a non-maximal halo, since in 
that case, the change of variable from z to ^ is not determined by the potential of Eq . (|3~3"1) 
but by an axisymmetric analogous of Eq.(|17|). Therefore the analytic expression given by 
Eq.(|36"D can be used only for maximal halos, i.e. po=p™ ax ■ 

The second family of distribution functions is a generalization of the axisymmetric loga- 
rithmic potential to the case of an asymptotically non-flat rotational curve, while preserving 



the property of Eq. (|35|). This is obtained for the axisymmetric power-law potential (16 
(models C3 and C4): 

w = ( ^ 0) - (37) 

The corresponding matter density is: 

n (P A- ^« R " (^ 2 + m 2 c + (l-(3q 2 )R 2 + [2-q- 2 (l + (3)]^ ^ 



Evaluating Eq. (|38|) for R = Rq, z = the parameter \l/ a can be expressed in terms of the 
density p . Eqs.(|37D and (^) are the axisymmetric generalization of Eqs. (^4j) and ( |23| ) 
and they posses the same properties already discussed in Section [111 A| . As for the case of 
the logarithmic potential, the rotational curve is obtained from Eq.(^) by substituting the 
radial coordinate r in Eq.(|25|) with the radial coordinate in the galactic plane R, and it is 
asymptotically falling with R if (3 >0 (model C3) and rising if j3 <0 (model C4). 

In analogy with the logarithmic case, also for the power-law model an analytic solution 
for the DF can be worked out |16| with the form of Eq. (^). This applies again only for a 
maximal halo. The analytic formulas for F + can be found for completeness in Appendix [A|. 
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1. Co-rotation and counter-rotation of the halo 



As mentioned before, the DF for an axisymmetric model is known up to an arbitrary 
odd component F— The DF we summarized above for the Evan's models all refer to the 
pure even component: they all have F_(e,L z ) = and posses no bulk rotation. The case 
F^(e,Lz) 7^ corresponds to the case of a rotating halo, where the number of particles 
moving clockwise around the axis of symmetry is different from that in the opposite sense. 

A family of DF's with bulk rotation can be studied by constructing an explicit example 
for F-. This can be done, starting from a generic F + , by considering the linear combination 



25, 26, 27 



where 



F-[c, L z ) — F right (e, L z ) — F ie f t (e, L z ) 



(39) 



right 



e,L, 



F + (e,L z ), v (j> >0 
0, < 



(40) 



and 



left 



(41) 



0, v<p > 0, 

F + (e,L z ), < 0, 

The distributions F right and Fi e f t describe the configurations with maximal with the 



same density profile as F + [ 25| . A DF with an intermediate value of can be obtained as 
a linear combination of F + and F_, or, equivalently, of Fi e f t and F r i g ht- 



F(e, L z ) = r]F right (e, L z ) + (1 - r))Fi eft (e, L z 



(42) 



The parameter rj ranges from 1 (maximal co-rotation) to (maximal counter-rotation) and 
is related to the dimensionless spin parameter A of the Galaxy by: A = 0.36 1?7 — 0.5 1 p6[ . 
In order to be consistent with the available extensive numerical work on galaxy formation, 
A should not exceed the value 0.05 implying 0.36 < rj < 0.64. For all the Evans 
model discussed in this Section, we will also study the co- and counter- rotating situations, 
adopting the two values r] = 0.36 and r] = 0.64 



D. Triaxial models 



The last class of models we wish to discuss is represented by the triaxial potential discussed 

1 



inRef. 0: 

T / x i 2, ( i y 2 z2 \ 

y {x,y,z) = --v ln\x + — + —J, 

which, for a maximal halo, corresponds to the DM density: 

vl Ax 2 + By 2 /p 2 + Cz 2 /q 2 



Pdm(x,V,z) 



4ttG (x 2 + y 2 /p 2 + z 2 /q 2 Y 



(43) 



(44) 



12 



where A = p 2 + q 2 — 1, B = 1 + q 2 — p 2 and C = 1+p 2 — q 2 . In Ref. Q the velocity 



DF f(v) of the system is approximated by a triaxial Gaussian with semiaxes equal to the 
velocity dispersions as obtained by the solutions of the Jeans equations: 



-2 = ^0 / 45 N 

(2 + 5)(p- 2 + q- 2 -l) 1 } 

vl = V ^ 2 - l) (46) 

+ 2(p~ 2 + q- 2 - 1) K ' 

« = (47) 



2(p- 2 + q 

when the Earth's position is on the major axis of the equipotential ellipsoid (models Dl and 
D2), and: 

-2 = V$P^ (4g) 

(2 + 5)(l + q- 2 -p- 2 ) { ] 

vl ~ V ° 2 ^" 2 -P" a ) (49) 

v 2 - V l¥-P-) (50) 



when the Earth's position is on the intermediate axis (models D3 and D4). In Eqs.( 
the quantity 5 is a free parameter that in the spherical limit (p = q = 1) quantifies the 
degree of anisotropy of the velocity dispersion tensor: 



IV. CONSTRAINING THE MODELS 

Once a given model is chosen for the velocity distribution function of the dark matter 
particles, the parameters of the model have to be fixed using observational data. Unfor- 
tunately, due to its "darkness" , all our knowledge of the halo is of indirect nature P3L 129" 



and it includes requirements on the circular rotational speed (constraints on its flatness and 
its value at the solar circle and in the outer regions of the Galaxy) as well as observational 
constraints on the local surface density of the disk and on the dispersion velocity of the 
bulge. In general, one should construct a composite model of the Galaxy where the DM is 
coupled to other components like the disk and the bulge, and the parameters describing the 
various components are varied independently requiring that the observational constraints 
are satisfied. 

We have already noticed that WIMP direct detection rates are particularly sensitive to the 
value of the rotational velocity v and the local DM density p (both evaluated at the solar 
circle). The procedure we follow in order to determine the allowed ranges for v and p i n 
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each galactic model is explained in the following. First of all, the experimental informations 
we use are: the allowed range for the local rotational velocity, the amount of flatness of 
the rotational curve of our Galaxy and the maximal amount of non-halo components in the 
Galaxy. The first information directly fixes the allowed interval for vq, irrespective of the 
galactic halo model. The other constraints are used in order to determine the allowed ranges 
for p , for each halo model. 
The allowed interval for vq is: 

w = (220 ±50) km sec" 1 (90% C.L.), (52) 



which conservatively relies on purely dynamical observations [30]. Proper motion mea- 



surements of nearby stars lead to similar estimates for the central value of Vq, with a 



significantly smaller uncertainty. However they are based on the assumption of circular orbit 
of the observed objects. For definiteness, we will use in the following three representative 
values for v , which correspond to its central value and to the boundaries of its allowed 90% 
C.L. range of Eq.(H): v = 170, 220, 270 km sec" 1 . 

For the three representative values of Vq we then determine the corresponding allowed 
ranges for p . For each halo model and for each value of Vq, we calculate, as a function of 
Po, two quantities: (i) the total amount of mass M vis in components other than the halo 
(e.g. : disk, bulge) which is necessary in order to match the given value of local rotational 
velocity Vq; (ii) the value of the rotational curve at a distance of 100 kpc from the center 
of the Galaxy: = v rot (R = 100 kpc). These two quantities are somewhat constrained 
from observations, even though their constraints are often obtained by using some degree of 
galactic modeling. We conservatively quote the following ranges |23|, f29"f : 

1 • 1O 1O M < M ms < 6 ■ 10 10 M o (53) 
0-8 -v << t °< 1.2 -v . (54) 

where M & denotes the solar mass. The first constraint limits the amount of non-halo com- 
ponents, while the second is a constraints of "essentially flatness" rotational curve: only 
galactic halo models which provide a rotational curve which does not deviate from a flat one 
more than 20% at 100 kpc are accepted. The constraint of Eq. (|54]) is compatible with the 
estimates of the galactic mass at large radii as obtained by the dynamics of satellites of the 
Galaxy [p^fl . 



The behaviour of M vis and v},®® as a function of po is shown if Fig. 1 (which refers to 
Vq = 220 km s _1 ), Fig. 2 (v = 170 km s _1 ) and Fig. 3 (vq = 270 km s" 1 ). Each line refers 
to a different halo model. The upper panels show that M vis is a decreasing function of p , 
since increasing the amount of dark matter in the Galaxy implies that less matter in other 
galactic components is required to support the rotational curve. On the other hand, the 
value of the rotational velocity in the outer Galaxy is totally supported by the dark halo, 
and it is larger for more massive haloes. 
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When the constraints expressed in Eqs. (|53[) and (gj) are simultaneously applied, an 
allowed interval for p may be derived for each halo model. 

The procedure outlined above may be used as a simple recipe for identifying the intervals 
for the local density parameter. However, some caution must be taken in the application of 
the bounds on M vis in Eq.(^). As already discussed in the previous Section, for the models 
of class A and class B we solve the Eddington equation in order to determine the velocity 
distribution function of dark matter particles. In this case we can take into account also the 
situation in which the local rotational velocity is only partially supported by the halo, and 
therefore we can apply the limits to M vis given by Eq. (|53D . Instead, the analytic models of 
class C and class D may be applied only to the extreme case of a fully maximal-halo, since 
the analytic formulae for f(v) are derived under the assumption that only the halo matter 
density is present. For these models (class C and class D) the only case we can deal with is 
that of a fully maximal halo, which corresponds to M v i s = 0. This represents a conservative 
upper limit for p . From Fig. 1-3 we can see that upper bounds to po obtained from the 
M vis = limit are only a few percent larger than what is obtained by imposing the lower 
limit of Eq. ([53"D : M vis = 1 • 1O 1O M . For consistency and simplicity, we will use M vis ~ as 
a lower limit also for models of class A and B. 

The allowed intervals for po that we obtain by imposing the bounds on M„ s and are 
listed in Table |T|. As discussed above, both values p™ m and p™ ax will be used in the next 
Sections to perform the modulation analysis of the DAMA/Nal experimental data, while 
only p™ ax will be used for the analysis of the data for models of class C and D. We stress 
that the reason for this stands in the fact that for the analytic models of classes C and D, 
we know f(v) only for a maximal-halo. 



V. RESULTS AND DISCUSSION 

In this Section we make use of the halo models described in Section III, with the choice 
of parameters shown in Table [H], to analyze the annual-modulation signal present in the 
DAMA/Nal data |2| in terms of relic WIMPs with purely coherent interactions. The pro- 
cedure is the one outlined in Section II. The results are presented as 3a annual-modulation 
regions shown in the plane £crg™ 1 ^ 1 r eon ' ) versus mw- All figures are divided in three panels 
which correspond to the cases: v = 170, 220, 270 km sec -1 . 

A general feature of all the models is that by raising the parameter t> the modulation 
region moves from the upper-right to the lower-left of the mw-£&w&ilr °^ plane. This is 
easily understood since, for a given DM density profile, higher values of vq imply higher 
values of p (through Eq. (|2])) and of the velocity ellipsoid =< ViVj > (through the Jeans 
equations); the experimental value of the signal and the measured WIMP-nucleus recoil 
energy are fixed by the data, therefore the modulation region moves downward because: 

dRdet (nucleoli) 



dE 



^ PO X ^scalar I 00 ) 
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(see Eq. (|J)) and moves to lower masses because: 

E R oc m w x < v 2 > . (56) 

On the other hand, when one compares the different density profiles that we have dis- 
cussed in the previous Sections, it is worth noticing that, for a given value of v , a stronger 
singularity in the galactic center lowers the value of p™ ax (in order to keep constant the mass 
integral) with the consequence that the lower part of the modulation region rises. So the 
smallest values of i&^i^ ^ are reached by the models with a less singular density profile. 

Figures 4-18 show the result of the analysis for the models with a spherically symmetric 
density profile (models AO-7, Bl-7). Models AO-7 have an isotropic velocity dispersion, 
while in models Bl-7 a degree of anisotropy in the velocity dispersion is introduced through 



the Osipkov-Merrit term of Eq . fl27|) . The amount of anisotropy is controlled by the (3q 
parameter of Eq. (^) and it has been fixed at the value: (3q = 0.4. This value corresponds 
to a radial anisotropy. In order to account for the uncertainty in the DM local density po, for 
each case two regions are given, the higher (lower) corresponding to po — P™ m (Po 10 ^), where 
p™ m and p™ ax have been obtained as discussed in Section IV and are reported in Table |TJ. 
The numerical values of the parameters summarized in the third column of Table |I| have 
been chosen in order to ensure the compatibility of the curves of Figs. [I], [| |3| with the 
constraints of Eqs.( [5T5||>3| ), discussed in Section [TV]. The Jaffe models A4, B4, which goes 
into the isothermal sphere when R c — > oo, have been calculated for the smallest allowed 
value of the core radius R c in order to examine the case of maximal departure from the 



usual scenario. As a consequence of this, for this model p™ m ~ p™ ax and in Figs. |8|,[X5| the 
upper and lower modulation regions are superimposed. 

The effect of radial anisotropy in the velocity dispersion tensor (v r > vg — v^), which 
occurs for the models of class B, may be seen by comparing Figs. 0-[Tl| with the corresponding 



Figs. |T2|-|T%|. As a general feature, a reduction of the modulation effect is expected, since 



the WIMPs phase space is depopulated along the direction of the Sun's velocity. This is 
confirmed by the fact that in most cases the modulation regions move upwards and widens, 
although the size of the effect can be small. The effect of radial anisotropy on the WIMP 
mass is more involved. In particular, the modulation regions for models Bl, B3, B5, B6 
extend to heavier WIMP masses compared to the corresponding isotropic cases, while for 
models B2, B4, B7 the region moves to smaller WIMP masses. 

As already pointed out, for the models belonging to classes C and D only the regions for 
Po=p™ ax are shown. As far as the axisymmetric models of class C are concerned, they are 



shown in Fig. [19], [22], ^g. In each model the flatness parameter q has been chosen in order 
to have the maximal flatness compatible with observations and with the positivity of the DF. 
The main effect induced by flatness is through the increase in the local density p , as can be 
seen in Table [TTT]. As a consequence, the modulation regions for these models reach values 
of ^sc^iar 01 ^ significantly below the levels of the spherical cases. In Figs. |2(| |23|, |26|, |29| the 



same models are shown with a co-rotation effect of the halo implemented through Eq.([42] 



16 



with 7] = 0.64, while in Figs. |21], 0, |27], ^ for the same models a counter-rotation effect 
with T) = 0.36 has been introduced. The main consequence of halo co-rotation is a decrease 
of the relative velocity between WIMPs and the Earth. The energy of WIMPs is therefore 
smaller and in order to produce the same recoil energy in the detector the WIMPs have to 
be heavier. This can be verified in the figures, where the modulation region of co-rotating 
models may reach very high WIMP masses, even higher than 200 GeV. By the same token, 
in counter-rotating models the modulation region is shifted toward lower masses. 

The peculiar shape of the modulation region of Fig. |20| deserves some comments. The 
two disconnected closed contours, which arise at different mw values, are indicative of the 
superposition in the WIMP phase space of two components with well separated r.m.s. ve- 
locities. This is exemplified in Fig. |35| for the models Bl (q = 1) and CI (q = l/y/2), where 
the contour plots of the corresponding DF's are plotted in the v-v^ plane (in the galactic 
rest frame) and in the w-w^ plane (in the Earth's rest frame). It is evident from the figure 
that, in flattened models, small \L Z \ orbits are depopulated compared to the spherical case, 
leading to two well separated populations with v </> > and < 0. These two components 
have the same temperature in the galactic rest frame, but develop different r.m.s. velocities 
when boosted in the Earth's rest frame. As shown in Fig. |35|, the relative weight of the 
two populations can be tilted towards slow or fast WIMPs in the two cases of a co-rotating 
or counter-rotating halo, respectively. This explains why the co-rotating model of Fig. [20 
develops a second minimum at high WIMP masses, and this is not observed in the cor- 
responding non-rotating or counter-rotating cases. We have numerically verified that the 
peculiar disconnected region at high WIMP masses reduces in size when the flatness param- 
eter is increased (since in this way the velocity distribution function becomes more similar 
to the non-flattened one), while it shifts towards lower masses when the core radius R c is 
decreased. 

We conclude the discussion of our results with the triaxial models shown in Figs. [31], 
32] , [33| , p4| where, to be definite, the same choice of parameters of Ref . |L7|] is adopted. For 
these models a general solution for the DF is not available. Only the velocity ellipsoid of Eq. 
(|4"5|-|4"8D is known, and it is used to fix the second moments of a non-isotropic Maxwellian. 
This explains why the shape of the modulation regions is quite similar to the standard case. 
In models Dl and D2 the Earth is assumed to be located on the major axis of the density 
ellipsoid, while in models D3 and D4 it is placed on the intermediate axis. Since in the two 
cases the Sun's position Rq is the same, po is higher for models Dl, D2 than for models D3, 
D4 (see Table |Hl[) . As a consequence, the modulation regions of Fig. |3l], [32| reach smaller 
values of ^(jg^ 1 r eon ' ) compared to those of Figs. |33], |4]. Models Dl and D3 (D2 and D4) have 
5 = —1.78 (5 = 16), which implies a radial (tangential) anisotropy of the velocity ellipsoid 
(see Eq. (pi]) ). Solving the Jeans equation for the potential of Eq.([|3|) in the spherical limit 
q = p = 1 (which corresponds to a non-cored isothermal sphere) leads to the relation JT]J : 

vl + vl = \v 2 . (57) 
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Note that this property is a consequence of the flatness of the rotational curve, and would 
not be true, for instance, including a core radius in the potential. Eq.floTD implies that, 
in this model, the tangential components of the velocity ellipsoid are fixed by Vq. As a 
consequence of this, radial anisotropy (v r > Vg — v^) corresponds to faster WIMPs and 
tangential anisotropy to slower WIMPs. This shifts the modulation regions towards smaller 
values of m w in Figs, 1ST], BBI and higher values of m w in Figs. This effect is 



sizeable in the case of tangential anisotropy, where the modulation region may extend up to 
m w ~ 270 GeV. 



The results of this section are summarized in Fig. 36, where all the modulation regions 
previously discussed have been plotted jointly. A convolution of all the regions may be 
indicative of the uncertainties in the determination of the WIMP modulation signal due to 
the modeling of the WIMP DF. As a final result, we show such a convolution in Fig. |37| , 
where a single curve in the plane mw^o'^^ * s plotted by collecting the information 
contained in the analyses of all the non-rotating models considered in this paper. The 
region is compared with the original annual modulation region obtained in Ref. @ for an 
isothermal sphere model of the galactic halo with rotational velocity v$ = 220 km s _1 and 
local dark matter density po — 0.3 GeV cm -3 . From Fig.[37| we see that the DAMA/Nal 
annual modulation result is compatible with WIMPs masses up to mw — 270 GeV and 
WIMP-nucleon cross sections in the interval: 10~ 10 nbarn < ^cig™ 1 ^ 1 r eon ' 1 < 6 x 10~ 8 nbarn 
when the uncertainties in the WIMP velocity DF is taken into account. Co-rotating models 
with maximal corotation can extend the mass range even further, up to mw ~ 500-900 
GeV, for cross section of the order few x 10~ 9 nbarn < £&^iar° n ^ £ 2 x 10~ 8 nbarn, as it 
can be seen, for instance, in Fig. EM 



VI. CONCLUSIONS 

In the present paper we have extended a previuos analyses of the DAMA modulation 
experiment for the case of a WIMP with a purely spin-independent coupling, by discussing 
in detail the implications on the results of the uncertainties on the dark matter galactic 
velocity distribution. We have studied a large number of viable models which deviate from 
the standard isothermal sphere in the matter density profile, in the presence of anisotropies 
of the velocity dispersion tensor and in effects of rotation of the galactic halo. The different 
models have been classified according to the symmetry properties of their matter density 
profile (or gravitational potential) and of the velocity distribution function. We have specifi- 
cally considered: (a) spherically symmetric matter density with isotropic velocity dispersion; 
(b) spherically symmetric matter density with non-isotropic velocity dispersion; (c) axisym- 
metric models; (d) triaxial models. 

The different models have then been used to re-analyze the DAMA/Nal 0-4 data collected 
by the DAMA/Nal Collaboration p[; in particular a total exposure of 57986 kg day, which 
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corresponds to 4 annual cycles, has led to the observation of an annual modulation effect. 
The hypothesis of WIMP annual modulation, already favoured in the previous studies |fj, 
H by using an isothermal sphere, is confirmed in all the investigated scenarios, and the 
effects of the different halo models on the determination of the allowed maximum-likelihood 
region in the WIMP mass and WIMP-nucleon cross-section have been derived. We can 
summarize that the DAMA/Nal annual modulation result is compatible with WIMPs masses 
up to my/ ~ 270 GeV and WIMP-nucleon cross sections in the interval: 10~ 10 nbarn < 
£ a scaiar !0n ' ) ^ 6 x 10~ 8 nbarn, when the uncertainties in the WIMP velocity DF is taken 
into account. When also co-rotation of the galactic halo is considered, the mass range 
extends further to my/ ~ 500-900 GeV, for cross section of the order few x 10~ 9 nbarn < 
^scalar 6011 " 1 ~ 2 x 10 -8 nbarn. These intervals quantify the extent of the annual modulation 
region for WIMPs with purely spin-independent couplings, as due to uncertainties in the 
phase space distribution function of galactic WIMPs. 
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APPENDIX A: AXISYMMETRIC MODELS 



In this appendix we give, for completeness, the analytic formulas of the DF's for the 
axisymmetric potentials of Eqs. ([33|,[37|) adapted from Refs. 
written in the reference frame of the Galaxy. 



15 , 16[ . All expressions are 



1. Logarithmic potential 

The DF for the logarithmic potential of Eq. fl33|) can be written as: 

F(e } L 2 z ) = (A + B)exp(^f]+Cexp(^] , (Al) 



where: 



o 



o x 1/2 p2 p2 i 

B = F 4*[±) j-£ *° 1 (A3) 



Rl + Rl 



C = F ,-^ J l J q \, \ (A4) 



and F = 0.47 GeV cm" 3 /^. 
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2. Power law potential 



The DF for the power-law potential of Eq. (|37|) for (3 > can be written as: 

F(e, L 2 ) = Ae 4 ^" 3 / 2 + Bi 4 ^^ 2 + Ce 2 ^" 1 / 2 , (A5) 

where: 



1„,2 



1 /v 



2 Vu! 



(A6) 



with C = R c / \j R 2 C + -Ro, while the velocity v\ = y\^? a \ is fixed through Eq. fl38"D , and with 
suitable normalizations, can be cast in the form: 



ui = 220 km sec 1 



Po 



Ra 



g( J R 2 + i? 2 )^+ 4 )/ 4 



v 0.47GeVcm-V 8.5 kpc J\j3\R Rp/ 2 R2(l + 2q 2 ) + i2g(l - /?g 2 ) ' 

(A7) 

The quantities ^4, S and C may be written as: 



A = F x 
B = Fx 
C = Fx 



R 



r ( 2 + 4 ^ a (0 + 2 )(!-l) 

220 km sec-V \R C J 2 3 /% 3 / 2 r(4//3 - 1/2)^ A g 2 J ' 



i? 



r(2 + 4//3) 



220 km sec- 1 ; Vi? c 7 2 3 / 2 vr 3 / 2 r(l/2 + + 2 \ 2 



vx 



2 ^n\ 2 r(2 + 2//3) 



220 km sec" 1 / Vi? c / 2 3 / 2 vr 3 / 2 r(l/2 + 2//?) 







2 - 



1 + /?' 



with Fx = 0.47 GeV cm _3 /ff. 

For p < Eq. (|A5|) still applies, with the following modifications (now < 0): 



— e 



where: 



A = Fx 



B = Fx 



C = Fx 



220 km sec" 



vx 



220 km sec" 



^(3/2-4//?) J_. 
i?J 2 3 /% 3 / 2 r(-l-4//?) PlP+ n g 2 ' 
-Re 



.220 km sec" 1 
with again Fx = 0.47 GeV cm~ 3 /vf 



r(l/2-4//3) 
i? c ; 2 3 / 2 7r 3 / 2 r(-l -4//3) 
^o\ 2 r(l/2-2//3) 
Rj 2 3 / 2 vr 3 / 2 r(-l -2//?) 



/3(/3 + 2)— . 



/3 



1 + 



- 2 



(A8) 
(A9) 
(A10) 



(All) 

(A12) 
(A13) 
(A14) 
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TABLE I: Values of the parameters for the spherically symmetric density profile of Eq.(|26|). 





a (3 7 a (kpc) 


NFW 
Moore 
Kravts 
Jaffe [ 


et 
o\ 
L4 


al. |19| 
et al. [20] 


13 1 20 
1.5 3 1.5 28 

2 3 0.4 10 

14 2 160 
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TABLE II: Summary of the galactic halo models discussed in Section III. The label introduced 



in the first column is used throughout the text to indicate each model in a unique way. For all 
the models the numerical values of the parameters which have been used in the calculations are 
given in the third column. The last column contains references to the models in the text. Models 
of class C have been analyzed also including co-rotation and counter-rotation of the halo through 
Eq. (H 



Class A: Spherical pdm? isotropic velocity dispersion 



AO 
Al 
A2 
A3 
A4 
A5 
A6 
A7 



Isothermal sphere 
Evans' logarithmic [15] 
Evans' power-law [16| 



Evans' power-law [16| 



Jaffe |14 
NFW HR 



Moore at al. |19] 
Kravtsov et al. ||2 



R c = 5 kpc 
R c = 16 kpc, = 0.7 
R c = 2 kpc, = -0.1 

Table | 

Table | 

Table | 

Table | 



Eq. 
Eq. 
Eq. 
Eq. 
Eq. 
Eq. 
Eq. 
Eq. 



Class B: Spherical pdm? non isotropic velocity dispersion (Osipkov Merrit, 0o = 0.4) 



Bl 

B2 
B3 
B4 
B5 
B6 
B7 



Evans' logarithmic 
Evans' power-law 
Evans' power-law 
Jaffe 
NFW 

Moore at al. 
Kravtsov et al. 



R c = 5 kpc 
R c = 16 kpc, = 0.7 
R c = 2 kpc, = -0.1 

Table | 

Table | 

Table | 

Table | 



Eqs. 
Eqs. 
Eqs. 
Eqs. 
Eqs. 
Eqs. 
Eqs. 



2|,||) 



Class C: Axisymmetric pdm 



CI 
C2 
C3 
C4 



Evans' logarithmic 
Evans' logarithmic 
Evans' power-law 
Evans' power-law 



Rc 
R c 



R c = 0,q= l/y/2 
R c = 5 kpc, q = l/y/2 
= 16 kpc, q = 0.95, = 
: 2 kpc, q = 1/V2, = 



0.9 
-0.1 



Eqs. 
Eqs. 
Eqs. 
Eqs. 



3§fMj) 



Class D: Triaxial p DM O (q = 0.8, p = 0.9) 



Dl 
D2 
D3 
D4 



Earth on major axis, radial anisotropy 
Earth on major axis, tangential anis. 
Earth on intermediate axis, radial anis. 
Earth on intermediate axis, tangential anis. 



5 = -1.78 
5 = 16 

S = -1.78 
5 = 16 



Eqs. 
Eqs. 
Eqs. 
Eqs. 



mm 

MM) 
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TABLE III: Allowed intervals of po obtained from the constraints on M v i s and v^t f° r t ne halo 
models summarized in Table O. The values of p™ ax and p™ 171 are used in the modulation analysis 
of the experimental data of Figs. |i|-fL8| for the models of class A and B, while only pQ ax is used for 



models of class C and D in Figs. 19-34. The value of p™ ax for the axisymmetric models of class C 
are not affected by the inclusion of a co-rotation or counter-rotation effect through Eq.(|42l). 





^0 


= 170 km sec 1 


^0 


= ZZu km sec 


^0 


= 270 km sec 1 


Model 


Po 


Po 


Po 


Po 


Po 


mas 

Po 


AO 


0.18 


0.28 


0.30 


r\ at 

0.47 


0.45 


0.71 


A 1 n i 

Al , Bl 


0.20 


0.42 


0.34 


0.71 


0.62 


1.07 


A O T"> O 

A2 , B2 


0.24 


0.53 


0.41 


0.89 


0.97 


1.33 


A3 , B3 


0.17 


0.35 


0.29 


0.59 


0.52 


0.88 


A4 , B4 


0.26 


0.27 


0.44 


0.45 


0.66 


0.67 


A5 , B5 


0.20 


0.44 


0.33 


0.74 


0.66 


1.11 


A6 , B6 


0.22 


0.39 


0.37 


0.65 


0.57 


0.98 


A7 , B7 


0.32 


0.54 


0.54 


0.91 


0.82 


1.37 


CI 


0.36 


0.56 


0.60 


0.94 


0.91 


1.42 


C2 


0.34 


0.67 


0.56 


1.11 


0.98 


1.68 


C3 


0.30 


0.66 


0.50 


1.10 


0.97 


1.66 


C4 


0.32 


0.65 


0.54 


1.09 


0.96 


1.64 


Dl , D2 


0.32 


0.50 


0.54 


0.84 


0.81 


1.27 


D3 , D4 


0.19 


0.30 


0.32 


0.51 


0.49 


0.76 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

p (GeV/cni 3 ) 



1 1.1 1.2 1.3 1.4 1.5 



FIG. 1: The quantities M v i s (upper panel) and v^f = v ro t(R = 100 kpc) (lower panel) calculated 
as a function of po f° r the halo models summarized in Table || and for fo = 220 km sec -1 . The 
different curves correspond, from left to right, to the following halo models: A4, AO, D3, A3, A6, 
Al, A5, Dl, A7, CI, C4, C3, A2, C2 (upper panel); A3, AO, D3, A5, Al, A6, A2, A4, C3, A7, 
C2, CI (lower panel: here C4 and Dl are not plotted because indistinguishable from A7). The 



horizontal lines indicate the constraints discussed in Section IV. The Bl-7 models have the same 
density distribution of the corresponding Al-7 models. D2 and D4 have the same distribution as 
Dl and D3, respectively. 



26 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 

p (GeV/cm- 3 ) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 

p (GeV/crrr 3 ) 

FIG. 2: The same as in Figure |l] for v$ = 170 km sec -1 (the ordering of all the different curves is 
maintained) . 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 

p (GeV/cm- 3 ) 



400 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 

p (GeV/crrr 3 ) 

FIG. 3: The same as in Figure |l| with vq = 270 km sec -1 (the ordering of all the different curves 
is maintained). 
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FIG. 4: Plot of the 3cx annual-modulation region in the plane ^^^° n ^ versus mjy using for the 
velocity distribution of WIMPs the isothermal sphere model (Model AO, see Table The three 
panels of the figure correspond to Vq=170, 220, 270 km sec -1 from left to right. Upper (lower) 



regions correspond to po=Po im (p™^) where p^ m and p™ ax are given in Table III 
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FIG. 5: The same as in Fig. || for the velocity distribution of model Al (see table [H]). 
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FIG. 6: The same as in Fig. for the velocity distribution of model A2. 
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FIG. 7: The same as in Fig. for the velocity distribution of model A3. 
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FIG. 8: The same as in Fig. [| for the velocity distribution of model A4. In this case p™ ax ~ />™ m , 
so upper and lower curves are not distinguishable. 



le-07 




le-10 



50 100 150 

m w (GeV) 



50 100 150 50 100 150 200 

m w (GeV) m w (GeV) 



FIG. 9: The same as in Fig. 0] for the velocity distribution of model A5. 



30 



le-07 




le-10 



50 100 150 50 100 150 50 100 150 200 

m w (GeV) m w (GeV) m w (GeV) 



FIG. 10: The same as in Fig. for the velocity distribution of model A6. 
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FIG. 11: The same as in Fig. |I| for the velocity distribution of model A7. 
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FIG. 12: The same as in Fig. with anisotropy of the velocity dispersion through the Osipkov- 
Merrit term and anisotropy parameter (3q = 0.4 (model Bl). The horizontal axis has been extended 
in the first panel. 
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FIG. 13: The same as in Fig. || with anisotropy of the velocity dispersion through the Osipkov- 
Merrit term and anisotropy parameter /3q = 0.4 (model B2). 
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FIG. 14: The same as in Fig. [7], with anisotropy of the velocity dispersion through the Osipkov- 
Merrit term and anisotropy parameter (5q = 0.4 (model B3). The horizontal axis has been extended 
in the first panel. 
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FIG. 15: The same as in Fig. |8|, with anisotropy of the velocity dispersion through the Osipkov- 
Merrit term and anisotropy parameter (3q = 0.4 (model B4). In this case pQ ax ~ Po' m , so upper 
and lower curves are not distinguishable. 
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FIG. 16: The same as in Fig. with anisotropy of the velocity dispersion through the Osipkov- 
Merrit term and anisotropy parameter (3q = 0.4 (model B5). The horizontal axis has been extended 
in the first panel. 
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FIG. 17: The same as in Fig. 10, with anisotropy of the velocity dispersion through the Osipkov- 
Merrit term and anisotropy parameter (3q = 0.4 (model B6). The horizontal axis has been extended 
in the first panel. 
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FIG. 18: The same as in Fig. 11, with anisotropy of the velocity dispersion through the Osipkov- 



Merrit term and anisotropy parameter /3q = 0.4 (model B7). 
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FIG. 19: The same as in Fig. || for the velocity distribution of model CI. Only the case po=p™ 
is shown. 
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FIG. 20: The same as in Fig. |19| including a co-rotation effect of the halo with r\ = 0.64. The 
horizontal axis has been extended in all panels. 
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FIG. 21: The same as in Fig. 19 including a counter-rotation effect of the halo with r] = 0.36. 
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FIG. 22: The same as in Fig. || for the velocity distribution of model C2. Only the case po=p™ 
is shown. 
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FIG. 23: The same as in Fig. |22| including a co-rotation effect of the halo with r\ = 0.64. The 
horizontal axis has been extended in the first panel. 
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FIG. 24: The same as in Fig. 22 including a counter-rotation effect of the halo with r] = 0.36. 
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FIG. 25: The same as in Fig. [| for the velocity distribution of model C3. Only the case Po=Pq l 
is shown. 
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FIG. 26: The same as in Fig. 25 including a co-rotation effect of the halo with rj = 0.64. 
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FIG. 27: The same as in Fig. 25 including a counter-rotation effect of the halo with rj = 0.36. 
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FIG. 28: The same as in Fig. [| for the velocity distribution of model C4. Only the case Po=Pq l 
is shown. 
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FIG. 29: The same as in Fig. 25 including a co-rotation effect of the halo with rj = 0.64. 
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FIG. 30: The same as in Fig. 28 including a counter-rotation effect of the halo with rj = 0.36. 
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FIG. 31: The same as in Fig. [| for the velocity distribution of model Dl. Only the case po=p™ 
is shown. 
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FIG. 32: The same as in Fig. |] for the velocity distribution of model D2. Only the case po=Po 
is shown. The horizontal axis has been extended in the first panel. 
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FIG. 33: The same as in Fig. |] for the velocity distribution of model D3. Only the case Po=Pq 1 
is shown. 
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FIG. 34: The same as in Fig. |] for the velocity distribution of model D4. Only the case po=Po 
is shown. The horizontal axis has been extended in the first panel. 
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FIG. 35: Contour plots of the DF's for models Bl and CI (see Table ||). From left to right, the 
different panels refer to models Bl (panels (a) and (e)), CI (panels (b) and (f)), co-rotating CI 
(panels (c) and (g)), counter-rotating CI (panels (d) and (h)). Upper panels are plotted in the 
Vff,-v plane, defined in the reference frame of the Galaxy, while lower panels are shown in the w^-w 
plane, defined in the reference frame of the Earth. Solid lines, big dashes, small dashes and dots 
correspond to growing values of the DF (in arbitrary units). The two disconnected closed contours 
which arise at different w values in panels (f), (g) and (h) signal the superposition in the WIMP 
phase space of two components with well separated r.m.s. velocities. 
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FIG. 36: Summary of the 3d annual-modulation regions in the plane £^scaiar° n ^ versus m\y, ob- 
tained by superimposing the results obtained with the velocity distributions of all the models 
described in Table O. For each of the models Al— 7 and Bl— 7 two regions are plotted, which refer 



to the two extreme values p™" and p^ ax shown in Table |TJ for the WIMP local density po- For 
models CI— 4 and Dl-4 only the regions which refer to po 



pmax are s h own _ 
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FIG. 37: 3a annual-modulation region in the plane £0g Calar versus mw, obtained by considering 
all the (non-rotating) galactic halo models discussed in this paper. The region is compared with 
the original annual modulation contour (shaded region) obtained in Ref. || for an isothermal 
sphere model of the galactic halo with rotational velocity vq = 220 km s _1 and local dark matter 
density po = 0.3 GeV cm -3 . 



41 



